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Abstract 

Decoupling via the Similarity Renormalization Group (SRG) of low-energy nuclear 
physics from high-energy details of the nucleon-nucleon interaction is examined for 
two-body observables and few-body binding energies. The universal nature of this 
decoupling is illustrated and errors from suppressing high-momentum modes above 
the decoupling scale are shown to be perturbatively small. 



1 Introduction 



The Similarity Renormalization Group (SRG) [T|[2|[3] provides a compelling 
method for evolving internucleon forces to softer forms [H[5] . While observ- 
ables are unchanged by the SRG's unitary transformations, the contributions 
from high-momentum intermediate states to low-energy observables is mod- 
ified by the running transformation. In particular, the SRG as implemented 
in Refs. [H[5] has the effect of partially diagonalizing the momentum-space 
potential to a width of order the evolution parameter. Because of this partial 
diagonalization, one anticipates a direct decoupling of low-energy observables 
from high-energy degrees of freedom. 

In Ref. [5], evidence for decoupling at low momentum was shown for the 
Argonne VI 8 [6] potential in calculations of phase shifts and the deuteron. 
In this letter, we extend the demonstration of decoupling to nucleon-nucleon 
(NN) potentials from chiral effective field theory (EFT) [7p8] and to few-body 
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nuclei up to A = 6 to verify its universal nature and to show quantitatively 
that the residual coupling is perturbative above the energy corresponding 
to the SRG evolution parameter. One might imagine that there is little to 
gain by evolving already-soft chiral EFT potentials to lower momentum, but 
recent work shows significant advantages for the nuclear few- and many-body 
problem [9]. This is consistent with the additional decoupling we find here. 

The practical test for decoupling is whether changing high-momentum matrix 
elements of the potential changes low-energy observables. Our strategy is to 
first evolve the initial potential Vnn with the SRG equations to obtain the SRG 
potential Vg, where s denotes the flow parameter of the transformation. Then 
we apply a parametrized regulator to cut off the high-momentum part of the 
evolved potential in a controlled way. This cutoff potential is used to calculate 
few-body observables and their relative errors. By varying the parameters of 
the regulator and correlating them with errors in the calculated observables, 
we have a diagnostic tool to quantitatively analyze the decoupling. 

In Section [21 we give background on decoupling in the SRG and the use of a 
smooth regulator to quantify decoupling. In Section [31 we analyze the relative 
error between phase shifts calculated with uncut and cut potentials and find a 
clean power law dependence at larger momenta that follows from perturbation 
theory. In Section [H we examine observables for the deuteron bound state 
and find the same perturbative behavior for the errors. Using the No-Core 
Shell Model (NCSM) pil|lllll2lll3|14j . we verify this behavior for few-body 
calculations in Section [5] and summarize in Section [61 



2 Decoupling in the SRG 

As in Ref. [1], we apply the similarity renormalization group (SRG) trans- 
formations to NN interactions based on the flow equation formalism of Weg- 
ner [2\. The evolution or flow of the Hamiltonian with a parameter s is a series 
of unitary transformations. 



where T^ei is the relative kinetic energy and H = Trei + is the initial Hamilto- 
nian in the center-of-mass system. Equation ([T|) defines the evolved potential 
Vg, with Trei taken to be independent of s. Then Hg evolves according to 



Hg = UgHUl = T,ei + Vg 



(1) 



s 



[Vs, Hg] , 



(2) 



ds 



with 



dUg 



(3) 



ds 
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Choosing rjs specifies the transformation, which is taken as the commutator 
of an operator, Gs, with the Hamiltonian, 

vis) = [G,,H,], (4) 

so that 

^ = [[G.,,H,],H,]. (5) 

Apphcations to nuclear physics to date in a partial-wave momentum basis 
have used Gs = Trei [1], but one could also use momentum-diagonal operators 
such as Tj,gj, Tj.gi, or the running diagonal Hamiltonian Hjj, as advocated by 
Wegner [2]. 

The source of decoupling is the partial diagonalization of the Hamiltonian by 
the SRG evolution. For the NN interaction, the flow equation ([5]) can be simply 
evaluated in the space of discretized relative momentum NN states [1]. For a 
given partial wave, with units where h'^/M = 1, we define diagonal matrix 
elements of momentum k as 



{k\Hs\k) = {k\Hn\k)=ek, (6) 

and 

(fc|T,ei|fc) = e, = . (7) 
With Wegner's choice, Gg = H^,, the flow equation for each matrix element is 

^^{k\Hs\k') = ^(e, + eu' - 2e,){k\Hs\q){q\Hs\k') . (8) 

By considering the trace of H^, one can show that off-diagonal matrix elements 
as a whole decrease in magnitude unless there are energy degeneracies [2] . 

If we take Gg = Tj-ei, the flow equation for each matrix element is 



j-^{k\Hg\k') = Y.{tk + efc' - 2t,){k\HM){q\Hs\k') 

= -{tk-tk'f{k\VsW) 

+ E(efc + ^k' - 2e,){k\Vs\q){q\Vs\k') . (9) 
q 

In this case, it can be shown that off-diagonal elements are not guaranteed to 
decrease monotonically if — e^' and — e^' have opposite signs (see Ref. [I5] 
for details). However, this does not happen in the range of s that has been 
considered in the nuclear case because of the dominance of the kinetic energy. 

This dominance can be used to find a semi-quantitative approximation for the 
flow of off-diagonal matrix elements by keeping only the first term on the right 
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side of Eq. ([9]) (or, equivalently, substituting {k\Hs\q) — (-k^kq)- Then the flow- 
equation apphed to individual off-diagonal matrix elements simplifles to 

^{k\H,\k') = ^{k\V,\k') ^ -{tk - tk^ f {k\Vs\k') , (10) 
as as 

which has the simple exponential solution 

{k\Vs\k') ^ (fc|K=o|A;')e-^(^'=-^^')' . (11) 

In Fig. [H we plot {k\Vs\k') ~ {k\Vs=o\k') and the approximation from Eq. ( fTTI) 
versus s for some representative off-diagonal points in two partial waves. In 
almost all cases the approximation gives a reasonable estimate of the mono- 
tonic decrease to zero; in the one exception there is a signiflcantly more rapid 
decrease than predicted. 
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Fig. 1. Absolute value of the matrix element {k\Vs\k') for a representative sampling 
of off-diagonal {k,k') pairs as a function of s, compared with the simple solutions 
from Eq. (jllj) . which are straight lines (they agree at s = 0). Two partial waves are 
shown. 

It is convenient to switch to the flow variable A = 1/s^^^, which has units of 
fm~^, since Eq. f lTTj) shows that A is a measure of the resulting diagonal width 
of Vs in momentum space. More precisely, the matrix {k\Vs\k') plotted as a 
function of kinetic energies fc^ and k'"^ will rapidly go to zero outside of a diag- 
onal band roughly of width A^, which is verifled by numerical calculations |16j . 
For momenta within A of the diagonal, the omitted quadratic part of the flow 
equation is, of course, essential, and drives the flow of physics information 
necessary to preserve unitarity. 

The tool we wih use to study decoupling is a smooth exponential regulator 
applied to the potential to cut off momenta above A: 

VUk, k') = e-('=^/^^)>,(fc, A;')e-(^'^/^^)" , (12) 
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where n takes on integer values. From the cut potentials we calculate observ- 
ables such as phase shifts and ground-state energies and compare to values 
calculated with the corresponding uncut potential. If there is decoupling be- 
tween matrix elements in a given potential (evolved or otherwise), we should 
be able to set those elements to zero in this systematic way and use the relative 
error in the observable as a metric of the degree of decoupling. By varying n 
we can identify quantitatively the residual coupling strength. 

The SRG transformations used here are truncated at the two-body level, which 
means that they are only approximately unitary for A > 3 and observables 
will vary with s. In these cases decouphng is tested by comparing cut to uncut 
potentials at a fixed s. All two-body observables calculated with the uncut 
Hg are independent of s to within numerical precision. The actual numerical 
error depends on the details of the discretization (e.g., the number and distri- 
bution of mesh points, usually gaussian) and on the accuracy and tolerances 
of the differential equation solver. While in practice we can make such errors 
very small, to avoid mixing up small errors we will also compare cut to un- 
cut potentials rather than to the unevolved (s = 0) potential for two-body 
observables. 



3 Phase Shift Errors 



In the upper-left panel of Fig. [2l we show results for the ^Sq phase shifts vs. 
energy calculated using the unevolved 500 MeV N^LO potential of Ref. [7] and 
the corresponding SRG potential evolved to A = 2.0 fm~^ and then cut using 
the regulator of Eq. (fT2|) with n = 8. We do not explicitly show results from 
uncut SRG potentials, because they are indistinguishable from the unevolved 
results. 

The qualitative pattern is that when the regulator parameter A is greater 
than A, there is good agreement of phase shifts from uncut and cut potentials 
at small energies and reasonable agreement up to the energy corresponding 
to the momentum of the cut, i^iab ~ 2ls? /m (with h = 1). When V^rg is 
cut below A, there is poor agreement everywhere and the phase shift is zero 
above this energy (e.g., above i^iab = 100 MeV for A = 1.1 fm~^). Thus the 
decoupling of high and low momentum means that we can explicitly cut out the 
high-momentum part of the evolved potential without significantly distorting 
low-energy phase shifts as long as we don't cut below A. Cutting out the 
high-momentum part of conventional nuclear potentials does cause distortions, 
which has led to the misconception that reproducing high-energy phase shifts 
is important for low-energy nuclear structure observables [S]. 

The quantitative systematics of SRG decoupling is documented in the other 
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Fig. 2. Phase shifts and relative errors in the So channel for SRG potentials evolved 
from the N^LO (500 MeV) potential of Ref. [7j. The upper-left graph shows the phase 
shifts vs. energy for the uncut A = 2fm~^ potential and several cut versions with 
n = 8. The other panels show the relative error as a function of the momentum cut 
parameter A at various energies E, A's, and n's, respectively. 

panels of Fig. [2], where we look at the relative error as a function of the cutting 
momentum using log-log plots. In these error plots, three main regions are 
evident. In the region below the A corresponding to the fixed energy, the 
predicted phase shift goes to zero since the potential has vanishing matrix 
elements, so that the relative error goes to one. Starting at A slightly above 
the value of A, there is a clear power-law decrease in the error. In between is 
a transition region without a definite pattern. 

We focus here on the power-law region. In the lower-left pane, we find that 
this decouphng starts with a shoulder at momenta slightly above A. This effect 
saturates when A becomes comparable to the underlying cutoff of the original 
potential. In the upper-right pane we see that the shoulder signaling the start 
of the power- law decrease is not affected by the energy, E. This holds for other 
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values of A and n. In the lower-right pane we vary the exponent of the regula- 
tor, n, which changes the smoothness of the regulator. The smoothness affects 
the slope of the power law and the fine details in the intermediate region, but 
does not change the position of the shoulder near A. As discussed below, the 
power-law behavior in the relative error signifies perturbative decoupling with 
a strength given by the sharpness of the regulator used to cut off the potential. 

We checked this decoupling behavior in different partial waves and for other 
N^LO potentials and found the same perturbative region in all cases. Repre- 
sentative partial waves are shown for various A values in Fig. [31 The potential 
in the S waves typically passes through zero for momenta in the region of 
A = 2fm~^ (see, for example, Ref. [H]), which might lead one to associate 
decoupling with this structure. The error plots for other partial waves that 
lack this structure show that it is a more general consequence of the SRG 
evolution. 




Momentum of Cut, A [fm '] Momentum of Cut, A [fm '] 



Fig. 3. The phase shift errors computed in select partial waves, using a regulator 
with n = 8. Other channels exhibit the same power-law dependence of the error for 
A > A. 

Indeed, the behavior of the errors in the decoupling region, where A > A, can 
be directly understood as a consequence of the partial diagonalization of the 
evolved potential. The calculation of the phase shift at a low-energy /c^ <C A^ 
will involve an integral over p of Vx^\{k,p). But the potential cuts off the 
integral at roughly ^ fc^ + A^ < A^, which means that we can expand the 
difference in the uncut and cut potentials: 

/ l.2n 2n \ 

5Vx,^{k,p) ^ V,ik,p) - V,Ak,p) ^ + X^^j V{k,p) . (13) 

Simple perturbation theory in 6V then predicts the dependence of the phase 
shift error to be 1/A^", which is the power-law dependence seen in Figs. |2]and 
[31 The accuracy of first-order perturbation theory is evidenced by the constant 
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slope of the error curves, which translates into perturbatively small residual 
coupling. 

The detailed dependence on the energy and A is not so trivially extracted. 
However, the weak dependence on i?iab < 100 MeV and strong dependence on 
A < 3fm~^ at fixed A seen in Fig. [2] implies that the integration picks up the 
scale A, so that the dominant error scales as (A/A)^". This is, in fact, observed 
numerically for intermediate values of A (e.g., for 1.8 fm~^ < A < 2.8 fm~^ 
when A = 3fm^^). 



4 Decoupling and Deuteron Observables 



To test the generality of the observations made for phase shifts, the same 
techniques were applied to other low-energy observables such as the deuteron 
binding energy, radius, and quadrupole moment. The binding energy and 
momentum-space wavefunction were computed using standard eigenvalue meth- 
ods. The computation of Qd and from the wavefunction uses |T7], 



/ dk 


V8( 


Jo 





and 



du{k) dw{k) 
dk dk 



+ ^kn^{k)'^] 
dk I 



(14) 



(15) 



where u{k) and w{k) correspond to the S and D components of the deuteron 
wavefunction respectively. We again computed relative errors in these observ- 
ables and, as shown in Fig. |4] for the energy and radius, the errors show the 
same behavior as observed for the phase shifts. That is, a power-law drop-off 
in the error begins at A just above A, with a slope determined by the sharpness 
of the regulator as given by n. The relative error for the quadrupole moment 
is not shown but is very similar to that for the radius. 

As with the phase shift, the analytic dependence of the error from cutting 
the potential can be estimated directly in perturbation theory. In this case, 
partial diagonalization of the potential means that the deuteron wave function 
has negligible momentum components starting slightly above A. This in turn 
validates the expansion in Eq. (fT3|l and the dependence of the errors on 1/A^". 
The numerical calculation of the error in perturbation theory is plotted in 
Fig. m and shows close agreement in the decoupling region A > A. 
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Fig. 4. The relative error vs. cut parameter A of the deuteron energy (left) and 
rms radius (right) with several values of the regulator parameter n indicated in 
the legends. In each case, the near-by solid line is the estimate of the error using 
first-order perturbation theory with Eq. 

5 Decoupling and Few-Body Energies with the NCSM 

The calculations described above have been only for two-particle systems. 
Using NCSM calculations of ground-state energies with the Many-Fermion 
Dynamics (MFD) code [18], we can test whether the high-energy decoupling 
behavior extends to few-body systems. In the present study, only NN interac- 
tions were considered, with the testing of decoupling with many-body forces 
deferred to a future investigation. However, we note that the general features 
of the SRG exhibited in Section [2] implies that off-diagonal matrix elements 
(with respect to energy) of the three-body force will be suppressed, with de- 
coupling as an expected consequence. 

We first verified that the decoupling behavior already observed using a direct 
calculation of the deuteron wavefunction is reproduced using the MFD. We 
then calculated a series of larger nuclei, including ^H, ^He, and ^Li, comparing 
results from uncut and a range of cut potentials evolved to different values of 
A. On the left panel of Fig. 0, the '^He ground-state energy is plotted versus the 
regulator parameter A for several different values of the SRG flow parameter 
A. Each of the plotted points is at a basis size A^max = 12; this is within several 
hundred keV of the energy from extrapolating to A^max = oo. A similar plot 
for ^Li is given in the left panel of Fig. using a basis size A^max = 8, which 
is within several hundred keV of the extrapolated energy for A = 1.5 fm~^ 
but still several MeV off for A = 3.0 fm~ . The uncut (A — > oo) energies vary 
for each A because the SRG evolution includes the NN interaction only; the 
closeness of the results for A = 2fm~^ and 3fm~^ for '^He is accidential (see 
Ref. [Hj for further discussion about the running of the energies). 
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Fig. 5. Calculations of the ^He ground-state energy using the NCSM. On the left 
is the energy obtained from the NCSM for potentials evolved to several different A 
values as a function of the cut (regulator) momentum A with n = 8. On the right 
is the relative error of the energy for the A = 2fm^^ case as a function of the cut 
momentum (with n = 8) for several different harmonic oscillator basis sizes. Also 
shown is the slope of the error in the decoupling region predicted from perturbation 
theory (dotted line). 
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Fig. 6. Calculations of the ^Li ground-state energy using the NCSM. On the left 
is the energy obtained from the NCSM for potentials evolved to several different A 
values as a function of the cut (regulator) momentum A with n = 8. On the right is 
the relative error of the energy for the same A's as a function of the cut momentum 
for the same A values but with two values of n. 

Both examples show that when the potential is cut with A comparable to A or 
lower, the converged energy is significantly different from the asymptotic uncut 
value, while it approaches that value rapidly as A moves above A. This means 
that, for smaller A, more high momentum matrix elements can be discarded 
without a loss of accuracy. This decoupling explains the greatly improved 
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convergence with basis size seen in the NCSM for corresponding A values [9]. 

The quantitative behavior of the relative error parallels that observed for two- 
body observables, as seen on the right panels of Figs. [5] and [61 In all cases, 
for a fixed value of A the power decrease in the error starting with A slightly 
above A is clearly seen, even though there are fewer digits of precision in the 
NCSM results (so the relative error is in the range 10~^-10~^ at best). The 
same perturbative residual coupling is seen for different basis sizes, with the 
slope given by the dependence although the onset of the decoupling 

region shifts to higher A until the calculation is near convergence (see Fig. 
Similar results are found for other nuclei and for other values of A. 

We repeated the calculations with other choices of the SRG generator, Gg 
including T^^i and Hd = Trd + Vd, where Vd is the (running) diagonal part of 
the bare potential. We found that these other choices for Gg do not alter the 
power-law behavior region of the previous error plots. This provides further 
evidence that the high- and low-energy decoupling results primarily from the 
partially diagonalized nature of the evolved potential. 



6 Summary 

The evolution of nucleon-nucleon potentials with the Similarity Renormal- 
ization Group decouples high-energy degrees of freedom from calculations of 
low-energy observables. By using a steep but smooth exponential function of 
the form exp[— (fc^/A^)"] with integer n to set interaction matrix elements with 
relative momenta above A smoothly to zero, we have found that the residual 
coupling follows a clear and universal behavior. Decoupling is achieved for A 
above A, the flow parameter for the SRG. The dependence on n was used to 
study the region of weak residual coupling, which was found to follow a power 
law predicted by leading-order perturbation theory. These results were shown 
to apply to NN phase shifts, several deuteron observables, and the ground- 
state energy of nuclei up to A = 6. Similar results are also found for other 
generators that are diagonal in momentum space. 

We emphasize that the regulator used here was only a tool for studying the 
effect of the SRG on the potential. No cuts of this nature have been used 
or are proposed for calculations with NN potentials. However, the most im- 
portant next step with the SRG is the evolution of three-body forces [T^ . 
which have significantly larger computational requirements. There is every in- 
dication that the SRG will induce the same decoupling for interactions that 
include three-body forces. These potentials will be projected on a momen- 
tum basis ordered by the kinetic energy of the states, and the SRG evolution 
equations will suppress matrix elements that are off-diagonal in that energy 
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basis [19] . Decoupling as observed here would allow a three-body computation 
to freeze irrelevant high-energy details during the evolution in a controlled 
way, which could significantly reduce the computational resources needed. 
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